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Gravitational-wave radiometry is a powerful tool by which weak signals with unknown signal 
morphologies are recovered through a process of cross correlation. Radiometry has been used, e.g., 
to search for persistent signals from known neutron stars such as Scorpius X-1. In this paper, we 
demonstrate how a more ambitious search—for persistent signals from unknown neutron stars— 
can be efficiently carried out using folded data, in which an entire ^year-long observing run is 
represented as a single sidereal day. The all-sky, narrowband radiometer search described here will 
provide a computationally tractable means to uncover gravitational-wave signals from unknown, 
nearby neutron stars in binary systems, which can have modulation depths of Ri0.1-2Hz. It will 
simultaneously provide a sensitive search algorithm for other persistent, narrowband signals from 
unexpected sources. 


I. INTRODUCTION 

Gravitational-wave radiometry is a technique by which 
two or more detectors are cross-correlated in order to 
identify signals in the form of excess coherence Hi- It 
is especially well-suited to situations where it is either 
impossible or impractical to carry out a matched filter 
search due to theoretical uncertainty and/or the vastness 
of the signal parameter space. LIGO @ and Virgo Q ra¬ 
diometer searches for persistent gravitational waves have 
yielded limits on the gravitational-wave strain from tar¬ 
gets including Scoroius X-I, the Galactic Genter, and 
Supernova 1987A 0. Radiometer searches are also sen¬ 
sitive to hot spots created from the superposition of many 
persistent sources [1-Q. By restricting the timescale 
of integration, the same method has been applied to 
search for relatively long-lived '^10-1000 s gravitational- 
wave transients [9|, e.g., associated with long gamma- 
ray bursts 0- Radiometry has also been proposed in 
the context of gravitational-wave astronomy with pulsar 
timing arrays |ll| . 

Due to its robustness and efficiency, there is strong 
motivation to extend the parameter space of radiometer 
searches. In this paper, we propose a method to extend 
the search for persistent narrowband signals to target 
all frequencies and all directions on the sky—not just 
ones associated with targets such as Scorpius X-1. An 
efficient all-sky search for persistent narrowband gravita¬ 
tional waves could facilitate the detection, e.g., of elec- 
tromagnetically quiet spinning neutron stars in binary 
systems. The challenge is to overcome computational 
challenges that make the search difficult. 

There is a significant burden from both data storage 
and computation time. In previous broadband radiometer 
searches for persistent gravitational waves Q, a Fisher 
matrix was used to characterize the covariance between 
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different patches of sky [I^. However, saving Fisher ma¬ 
trices for many frequency bins (and for many “jobs” as¬ 
sociated with stretches of science quality data) creates a 
burdensome data storage problem. The computational 
cost of analyzing a full year of data to look at thousands 
of frequency bins is both inefficient and potentially pro¬ 
hibitive in terms of computation time. 

Recent work [I^ shows that by “folding” data into a 
single sidereal day, it is possible to represent an entire 
year of radiometer data with just one twenty-four long 
spectrogram. We show that this several-hundred-fold re¬ 
duction in data volume facilitates efficient searches for 
persistent narrowband signals. 

The remainder of this paper is organized as follows. In 
Section m we discuss the motivation for narrowband all¬ 
sky searches. Section Hill describes how previous work on 
folded data [l^ and searches for f«day-long signals 0 
can be combined to design an efficient all-sky narrowband 
radiometer. In Section IlYl we demonstrate an all-sky 
narrowband radiometer search using folded data. Finally, 
in Section |Vl we discuss the astrophysical implications of 
these results. 


II. MOTIVATION 

One of the primary targets for an all-sky narrowband 
search for persistent gravitational waves is spinning neu¬ 
tron stars in binary systems. The motivation for such 
a search was recently laid out in [0. In a variety of 
scenarios, accretion from a companion star can induce a 
time-varying quadrupole moment, which, in some cases, 
may persist after accretion abates. For example, persis¬ 
tent localized mass accumulation may occur due to mag¬ 
netic fields [0 depending on unknown details of neu¬ 
tron star physics |17H^ . Magnetic fields may also in¬ 
duce deformations in the stellar interior (2^. Rotational 
instabilities such as r-modes may be sustained through 
accretion 0.0. 

When the neutron star in a binary is a pulsar, it is 
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possible to carry out an optimal search using matched 
filtering. However, if it is electromagnetically quiet, then 
the space of possible signals becomes prohibitively large 
for a fully coherent search. When at least the sky location 
is known, one can apply a targeted radiometer search [I|, 
or other semicoherent methods [l^ . [25l - l^ . If the sky 
location is unknown, the problem is significantly more 
challenging and the signal might be missed or vetoed as 
a noise artifact by non-specialized pipelines [l5l |. 

The search technique proposed here is highly comple¬ 
mentary to the semicoherent method described in |15l . 
[ 2 ^. The “TwoSpect” algorithm employs a canon¬ 

ical signal model in which a neutron star emits periodic 
waves, modulated by binary motion. The radiometer, 
on the other hand, employs only minimal assumptions; 
namely, that the signal is persistent and narrowband. 
By incorporating additional assumptions about the sig¬ 
nal model, a tuned search gains sensitivity. The advan¬ 
tage of the radiometer, in contrast, is that it is extremely 
robust since it relies only on excess coherence in two or 
more detectors. Indeed, while we focus on spinning neu¬ 
tron stars in binary systems, the narrowband radiometer 
is sensitive to any persistent narrowband source, whether 
or not it is modulated by binary motion. 


III. METHOD 
A. Folding 

We utilize the idea of folded data as described in 
which we reformulate in the language of recent radiome¬ 
ter developments [3, [ifl [3 . During an observing 

run, typically lasting a few months to a few years, data 
is collected from two spatially separated detectors. (This 
discussion straightforwardly extends to three or more de¬ 
tectors, but our presentation focuses on the two-detector 
case for simplicity.) We parse the data for which both de¬ 
tectors are simultaneously operational into sidereal days 
(23 hr, 56 min, 4 s). Every sidereal day of radiometer data 
can be represented using a complex-valued estimator M- 

( 1 ) 

aa, it; /) = i /)P'(t;/). (2) 

Here, (t; /) are spectrogram indices: t is the start time of 
each data segment and / is the frequency of the Fourier 
transformed data associated with that segment. Each 
s/(i;/) refers to the discrete Fourier transform of the 
strain data from detector I. The variable A/” is a Fourier 
normalization constant, a^it;/) is an estimator for the 
uncertainty associated with and each Pj{t;f) repre¬ 
sents the strain auto-power measured in detector I. (The 
prime denotes that the auto-power is calculated as an 
average of neighboring segments.) The direction of the 


source is encoded in the phase of 2)(t; /), which depends 
on the time delay between the two detectors [30l| . 

Since the signal-induced phase delay and detectors’ an¬ 
tenna factors are periodic on the timescale of a sidereal 
day, it is possible to define folded estimators by summing 
data over many sidereal days: 

k k 

( 3 ) 

■ ( 4 ) 

Here, k runs over all the sidereal days of the observing 
run. 

Following [ 33 , we define complex signal-to-noise ratio 
for folded data: 

The variable p^°*‘^(t; /) is a complex-valued spectrogram; 
see Fig. [TJ Gravitational-waves induce excess p(t;/), 
which appears visually as brighter than usual spectro¬ 
gram pixels. 

By folding a year of data into a sidereal day, the volume 
of data is reduced by a factor of ss365, which dramati¬ 
cally reduces both storage needs and computation time. 
And, if we restrict our attention to persistent narrowband 
signals, then the folding operation is lossless in the sense 
that there is no more information (useful in a search for 
persistent signals) in a year of cross-correlated data than 
there is in a folded day. 


B. Radiometry 

By specializing the methods from 0 to consider ap¬ 
proximately monochromatic signals lasting exactly one 
sidereal day, and then applying these methods to folded 
data 0 , it is possible to construct a very efficient algo¬ 
rithm for the detection of electromagnetically quiet bi¬ 
nary neutron stars (and other narrowband sources not 
matching a standard continuous wave model). Note that 
when we refer in this paper to signals as “approximately 
monochromatic,” we mean that they fit within single 
Ril Hz frequency bin. Such signals include neutron stars 
in binary systems, which are relatively broadband com¬ 
pared to isolated neutron stars. However, the signal is 
narrowband in the sense that it is contained within just 
one radiometer frequency bin. In reality, the modulation 
depth of neutron stars in binaries can range from ssO.l- 
2 Hz [ 1 ^, so a more nuanced definition of “narrowband,” 
employing variable frequency bin sizes, is appropriate. 
However, this is beyond our present scope. 

Radiometer searches can be cast as a pattern recog¬ 
nition problems in which an algorithm looks for statisti¬ 
cally significant clusters of pixels Q. As we are focused 
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here on approximately monochromatic signals, the opti¬ 
mal pattern recognition algorithm is a simple sum over 
the frequency bins of In the language of seedless 
clustering [^, , each signal template corresponds to a 

different row in the spectrogram. We assume that 
the signal starts at the beginning of the spectrogram and 
finishes at the end. To assume otherwise would imply 
an unphysical signal that starts and stops with a period 
matching Earth’s sidereal day. Likewise, it would not 
make sense to allow for any frequency evolution since 
this would imply that the signal repeats this evolution 
with the period of a sidereal day. 

The following expression for the detection statis¬ 
tic SNRtot(/|^) is specialized from 0 to focus on 
monochromatic signals: 


SNRtot(/|^) = 


Re 


g(2«/n.Ax(t)/c)pfold(^. 




1/2 


( 6 ) 

Here, 0 , is the direction of the source. Ax is the differ¬ 
ence in detector position, and c is the speed of light. 
The phase factor “points” the spectro¬ 

gram toward the source by rotating the phase angle of 
to zero so that the observed signal-to-noise ratio is 
real and positive 00. The variable ei 2 {t\fl) is a time- 
dependent efficiency factor characterizing the fraction of 
gravitational-wave power measured due to the non-unity 
antenna factors of gravitational-wave interferometers: 


A 


Here, Ff{t\Cl) is the antenna factor for detector I 
and H = -f, X are the different polarization states. As the 
Earth rotates, the detectors become more/less favorably 
aligned relative to the source. The sum over t is carried 
out over all the (typically ~60s) segments in a sidereal 
day. For additional information pertaining to Eq. |6l the 
interested reader is referred to [ij and reference therein. 

For the sake of compact notation, Eq. [5] assumes that 
detector noise is approximately constant over the span 
of a sidereal day. However, we note that it is straight¬ 
forward to generalize to the case of non-stationary noise 
by weighting each data segment appropriately using a 
standard inverse variance weighting method |34l | . 

Previous work [3, [sol - fs^ has shown that Eq. [5] and 
Eq. [S] can serve as the basis for “all-sky” searches for 
which the signal time and direction are unknown a priori. 
Many directions in the sky can be considered simultane¬ 
ously using parallelized computer code run on suitable 
processors, e.g., GPUs and multi-core CPUs 0. For 
each observation, we record SNR™“—the maximal value 
of SNRtot for all the directions: 


SNR“r(/)=max SNRtot(/|f^) 
n L 


( 8 ) 


Eq. [6] assumes a fairly constrained signal model: ap¬ 
proximately monochromatic signals, which start at the 


beginning of and hnish at the of a sidereal day. Thus, 
it is possible to search every s^ location (given some 
resolution) using a HEALPix grid. In this work, 
we consider a number of sky locations chosen so as to 
be spaced at angular separations less than the diffraction 
limited resolution of the detector network. For the LIGO 
network. 


c 1 /lOOOHzk 

/ |Af| V / J 


(9) 


In the demonstration described below, we use HEALPix 
to efficiently sample the sky with 3072 tiles. 


C. Comparison with targeted searches 

In order to highlight the usefulness of this framework, 
it is useful to compare it to the procedure used for previ¬ 
ous narrowband targeted radiometer searches 0 , 11 , ap¬ 
plied, e.g., to a single direction such as Scorpius X-1. In 
such searches, one analyzes data in RiSOO-lOOOO s-long 
“jobs,” defined as data stretches during which coinci¬ 
dent data is available for both detectors. During initial 
LIGO’s fifth science run, there were Ri20000 jobs. For 
each job, one calculates cross-power (Eq. [T|) and auto¬ 
powers (Eq. [2]) by optimally combining the results from 
each of many (typically) 60 s segments in each job. Then, 
during post-processing, the results from each job are op¬ 
timally combined to produce estimators analogous to 
and (Tp, but obtained by integrating over the entire ob¬ 
serving run. 

This procedure is repeated separately for each sky lo¬ 
cation. Using this scheme, each sky location requires the 
generation of tens of thousands of output files, and so 
there is significant overhead associated with the analysis 
of more than a few sky locations. Moreover, covariances 
between different patches on the sky—characterized by a 
Fisher matrix—lead to subtleties when interpreting the 
significance of an outlier 0| . In order to understand these 
covariances in this framework, one must therefore pro¬ 
duce also a Fisher matrix for every frequency bin. 

In the new paradigm we propose here, the data are 
first combined into a single sidereal day, which can be 
used to look in every sky direction. This data is small 
enough to be easily manipulated on a typical computer. 
Once the data are loaded, we can search over every sky 
location all at once without having to read or write addi¬ 
tional files, thereby eliminating a significant bottleneck. 
Also, this method eliminates the need to store Fisher ma¬ 
trices. This is because the covariances between different 
sky locations are automatically encoded into the factors 
of and e(t|H). 

Thus, one can use numerical simulations, in which p 
is drawn from a Gaussian distribution, in order to de¬ 
termine the false-alarm probability of obtaining a given 
value of SNR™“. (The assumption of Gaussianity is 
justified by invoking the central limit theorem, and it 
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has been born out in previous measurements H^-) 
An example search, analyzing 750 Hz of bandwidth with 
A/ = 1 Hz resolution (see Fig.[T]), and scanning over 3072 
Healpixels, can be carried out in < 100 s using a single 
eight-core CPU. This computation time does not include 
the modest time required for folding the data. 

IV. DEMONSTRATION 

Our goal now is to demonstrate the recovery of a per¬ 
sistent narrowband signal in folded data without a priori 
knowledge of the sky location or frequency of the sig¬ 
nal. To begin, we generate twenty days of simulated 
data for the LIGO Hanford and LIGO Livingston de¬ 
tectors operating at design sensitivity. The data con¬ 
sists of Gaussian noise plus a simulated monochromatic 
signal with frequency / = 600 Hz and strain amplitude 
ho = 1.5 X 10“^^, located at (ra,dec) = (18.5hr,39°). 
The signal is circularly polarized. For our purposes, a 
monochromatic signal is a reasonable proxy for a signal 
modulated by binary motion (discussed in Section |TT]) so 
long as the binary modulation is < 1 Hz. 

We generate coarse-grained spectrograms with a res¬ 
olution of (79s, 1 Hz) (with non-overlapping segments). 
This frequency bin width is relatively well-matched to 
the observed modulation depth of known neutron stars in 
binary systems, which can range from RiO.1-2 Hz [l^. (A 
more systematic study should employ variable bin widths 
to fully cover parameter space.) We study a band ranging 
from 500-1250 Hz, where many neutron stars in binary 
systems are expected to rotate [13 . Next, the data are 
folded into one sidereal day following the procedure out¬ 
lined in Section nni As the data are folded, we analyze 
the integrated spectrogram with every accumulated day 
in order to study how the signal grows with time. Fol¬ 
lowing the procedure described in Eqs. [6]and|8l we scan 
3072 Healpixels and record the maximum signal-to-noise 
ratio at each frequency SNR]^“(/). 

The results are summarized in Fig. [1] In Fig. [Iji, we 
show the real part of e(2'^®/sl'Aa:/c)pfoidj-^. represent¬ 
ing twenty days of data folded into a single sidereal day. 
Though the signal is virtually impossible to see with the 
naked eye in just one day of data, it is visible in this im¬ 
age at 600 Hz since the signal-to-noise ratio grows with 
continued integration. Fig.[TjD shows the recovered signal 
(SNRJ^“ Ri 12) using just one day of data. The excellent 
match in the time-dependent modulation suggest that 
the sky location is well matched. Indeed, the best fit di¬ 
rection (ra, dec) = (18.3 hr, 38.9°) is the closest possible 
pixel to the true source location. 

In Fig.[T]:, we show how the signal-to-noise ratio of the 
600 Hz signal (blue) grows as expected like the square 
root of the total observation time. The dashed red curve 
shows the next-loudest frequency bin, which exhibits no 
tendency to grow or decline with continued integration. 
In Fig. [T]i, we show the signal-to-noise ratio spectrum 
for the spectrogram in Fig. [TJi, maximized over search 


directions. That is, for each frequency bin, we scanned 
the entire sky and plot (in red) SNR”qj’° for the bright¬ 
est patch of sky as a function of frequency. The blue 
curves indicate typical one-sigma fluctuations for noise. 
The 600 Hz signal is visible as a dramatic spike. There 
is a slight tendency toward higher values of SNR with 
increasing frequency as angular resolution improves and 
there are more independent patches of sky. 

From Fig. [1] we show that it is possible carry out a 
computationally efficient search for persistent narrow- 
band signals from all directions in the sky and at all 
frequencies. Once the data have been folded, the entire 
search takes less than two minutes to carry out using a 
single eight-core GPU. 

By studying the background distribution of 
SNRj((j’'(/), we can estimate how the sensitivity of 
the radiometer search differs for a directed and an 
all-sky search. Naively, we expect the all-sky sensitivity 
to be worse, but only slightly, since the additional trial 
factors incurred by looking in many directions will only 
slightly increase the detection threshold given the ex¬ 
pected rapidly falling background distribution. In order 
to identify a detection candidate with false alarm proba¬ 
bility 1%, the targeted search requires SNR]^“(/) > 4.5 
while the all-sky search requires SNR"qj’°(/) > 5.8 for 
the bandwidth and spectral resolution considered here. 
This corresponds to a change in strain sensitivity of just 
(5.8/4.5)1/2 - 1 Ri 14%. 

Upper limits can be set for each frequency bin inde¬ 
pendently, and so there are fewer trial factors to apply 
compared to the detection calculation, which must ac¬ 
count for the existence of hundreds of frequency bins. 
On average, the targeted search will set limits corre¬ 
sponding to SNR]^“(/) Ri 1.7 whereas the all-sky search 
will set limits corresponding to SNRj((“(/) rs 4 (around 
/ R:! 1000 Hz where we scan r:!3072 sky positions). Thus, 
the all-sky limits are expected to be Ri50% higher in 
strain than limits obtained from a targeted search. These 
estimates support the expectation that the all-sky search 
is only slightly less sensitive than the targeted search. 


V. CONCLUSIONS 

Rotating neutron stars in binaries are a very promis¬ 
ing candidate for gravitational-wave detection by second- 
generation detectors like Advanced LIGO/Virgo. Previ¬ 
ous work (23 has enabled the exploration of these promis¬ 
ing sources, though, no detection candidates have been 
found with initial LIGO data [l3|. Our presentation of 
an all-sky, narrowband radiometer using folded data is 
complementary since it makes only minimal assumptions 
about the source. In this way, it should be sensitive not 
only to neutron stars in binary systems, but also to other 
sources, which emit approximately narrowband gravita¬ 
tional waves and which may or may not conform to the 
canonical model of an isolated neutron star. 

A back-of-the-envelope calculation provides promising 
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FIG. 1: Narrowband signals in folded data. Top-left: a {5t,Af) = (79 s, 1 Hz) spectrogram of the real part of 

g(27i-j/n.Aa:/c)pfoid^^. Consisting of twenty days of data, folded into one sidereal day. The data consist of Advanced LIGO 
Monte Garlo noise plus a simulated signal with amplitude ho = 1.5 x 10“^^ and frequency 600 Hz. The source is assumed 
to be face-on. The spectrogram is “pointed” in the correct direction using the appropriate phase factor to make the signal 
positive; see Eq. [T] The simulated signal is detectable in just one day of data, but we show the result of 20 days of folding 
to make the signal visible by eye. The top-right panel shows the recovered signal in just one day of data. The signal is easily 
identified without any prior assumption about the source location or frequency. The bottom-left panel shows how the folded 
SNR grows with the addition of folded data. Blue shows the result for the bin with the signal while red is the next loudest 
bin. The bottom-right panel shows the maximum SNR (scanning over the entire sky) as a function of frequency for one day of 
data (red). The injected signal is evident as a SNR « 12 spike at 600Hz. Typical fluctuations due to pure noise are indicated 
with blue. It is interesting to note that there is a slight trend toward higher SNR with increasing frequency due to increasing 
angular resolution. 


results. In Subsection IIVI we showed that a /iq = 1-5 x 
10“^^ signal at / = 600 Hz can be easily recovered in one 
day of data. Scaling to one year of data, this becomes 
/iQ ~ 2 X 10“^®. In terms of neutron-star ellipticity, this 
sensitivity can be restated [s^ as 


e 


1 X 10“® 



10^® gcm^ 
1 


/ r \ /600Hzy 

( 10 ) 


where /3 is an orientation factor, G is the gravitational 
constant, r is the distance to the source, and I is the 


moment of inertia. (One arrives at a similar rough esti¬ 
mate by scaling the upper limits from targeted radiome¬ 
ter analyses with initial LIGO data @ to Advanced LIGO 
sensitivity.) This is an order of magnitude above the ex¬ 
pected strain at / = 600 Hz from Scorpius X-1 expected 
from accretion-torque balance [2^ , suggesting that a sig¬ 
nal might be detectable from another somewhat closer 
neutron star—especially in light of the fact that addi¬ 
tional optimizations are possible, and only a hint of a 
signal is needed to trigger a coherent matched filtering 
search, which can boost the significance. 
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